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The isospin violating J — > 07t°t; decay is thought to be dominated by the mixing of the /o(980) 
and ao(980) scalar resonances. We make a theoretical evaluation of the J/^> — > (j>-K r) decay extending 
our own previous model for other J/'fy decays into one vector meson and two pseudoscalars using 
the techniques of the chiral unitary approach. The scalar resonances are dynamically generated 
through the final state interaction of the pseudoscalar mesons implementing unitarity from the 
lowest order ChPT amplitudes. Besides the direct J/ty(f)PP vertex, other mechanisms like the 
sequential exchange of vector and axial- vector mesons are shown to be important in order to obtain 
the actual strength of the mixing. We get a very good agreement with the Tv°ri invariant mass 
distribution and branching ratio with recent BESIII data. Quantification of the /o(980) — ao(980) 
mixing is done and compared with results from several experiments. 



I. INTRODUCTION 

The nature of the scalar mesons has been the subject of 
controversial debate for more than thirty years and a fi- 
nal consensus has not been reached so far. The problem 
stems from the fact that it is very difficult to accom- 
modate them within the oo picture [l[ . Inspired in the 
four quarks interpretation [2tJ] a second kind of molecule 
picture was developed, in which only mesonic degrees of 
freedom were considered Q. In the past decade new 
light has been shed by the chiral unitary coupled chan- 
nel approach [H-fl] which successfully describes meson- 
meson interaction in coupled channels up to energies far 
beyond the limit of applicability of standard chiral per- 
turbation theory (ChPT). With the only input of the 
lowest order ChPT Lagrangians, the implementation of 
unitarity in coupled channels and the exploitation of the 
analytic properties of the scattering amplitudes, many 
hadronic resonances appear dynamically as poles of the 
unitarized scattering amplitudes, e.g. [9l4l8|. These are 
usually called molecules or dynamica lly gen erated reso- 
nances, like the /o(980) and ao(980) |19h23| |. of interest 
for the present work. 

One interesting issue of the ao and /o resonances is that 
they can mix violating isospin, as it was first suggested 
in ref. 24]. This mixing is expected to help clarify the 
nature of these resonances. It is thought [HiHH to be en- 
hanced and only relevant in the narrow 8 MeV region be- 
tween the K + K~ and K°K° thresholds, (987.4 MeV and 
995.2 MeV respectively). Several reactions could be sen- 
sitive to the , fp — ao mixing, for example pn — > dr/Tr (26| , 
7P — > pTr°ri [27J or n~p — » n°rin [28[. Very recently 
the isospin violation in ?y(1405) decays into /o(980)-7r° 
have been experimentally measur ed [201] and theoretically 
studied in a simultaneous work 30] to the present one 
and using also techniques of the chiral unitary approach. 

Until very recently [3l[ there has been no conclusive 
experimental information on this mixing. In ref. [3l| 
the 7r°T7 invariant mass distribution of the J/*f? — > <fiir°ri 
decay has been measured experimentally by the BESIII 
collaboration in Beijing and a clear peak in the window 



between the K + K~ and K°K° thresholds has been ob- 
served. The branching ratio of the J/\& —> 0/o(98O) — > 
0ao(98O) — !• 4>ir°r) is determined in [31| to be of the or- 
der of 3.3 x 10 -6 with a significance of about 3.4ct which 
represents the most significant evidence of the isospin 
violating /o — ao mixing to date. From the theoretical 
side a first estimation of the mixing in the J/^f — > (f>ir°ri 
decay was done in ref. (32| using a simple model with ex- 
plicit scalar resonance prop agators coupling to the pseu- 
doscalar pairs. In ref. [33[ the J/\& — > 07r°77 was calcu- 
lated from the implementation of the final rescattering 
of the two pseudoscalars in the direct J/^ipPP vertex 
using the chiral unitary approach. 

In the work of ref. [34| we developed a very detailed 
and thorough model (one of the most comprehensive to 
date) for the decay of the J/^> into a vector meson and 
two pseudoscalars, in particular for the isospin conserv- 
ing decays into unr + ir~ , 4>tt + it~ , loK + K~ and <j)K + K~ . 
Besides the direct J/^<j)PP vertex and the implementa- 
tion of the final state interaction of the two pseudoscalars 
other new mechanisms were considered in ref. [34j where 
the J/vE" decays into a vector or axial- vector and a pseu- 
doscalar meson and the vector or axial-vector subse- 
quently decays into the final vector and another pseu- 
doscalar. We call these mechanisms sequential vector and 
axial- vector exchange. In this sequential mechanisms also 
the final pseudoscalar-pseudoscalar interaction is imple- 
mented from which the scalar resonances are naturally 
generated dynamically. We found a very good reproduc- 
tion of the experimental mass distribution also in the 
(t(500) and /o(980) regions, not trivial at all. 

In view of the power of the model of ref. [34j |. in the 
present work we are going to extend it to the J/\& — > 
rjm°r] decay in order to analyze the /o — ao mixing. The 
isospin violation will be implemented essentially from 
leading order quark mass differences and electromag- 
netism, similarly as in ref. (33^ . In practice it will come 
essentially from the difference between the masses of the 
charged and neutral kaons and pions in the loops and the 
meson-meson scattering amplitudes. One of the main dif- 
ferences with ref. (33| is the new sequential mechanisms 
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of ref. which, advancing some results, turn out to be 
important due to strong interferences with the dominant 
loops from direct production. 



II. THE f - a MIXING THROUGH 
MESON-MESON UNITARIZATION 

In the literature several unitarization procedures have 
been used to obtain a scattering matrix fulfilling exact 
unitarity in coupled channels, like the Inverse Amplitude 
Method [1 iH or the N/D method Q. In this latter 
work the equivalence with the Bethe-Salpeter equation 
used in [l9| was established. In the present work we use 
the Bethe-Salpeter equation which leads to the following 
unitarized amplitude in coupled channels 



t= [l-VG\- l V. 



(1) 



Diagrammatically Eq. ([T]) is equivalent to resum the 
series expressed by the thick dot in Fig. [2] In Eq. ([T]) G 
is a diagonal matrix with the I— th element, G;, being the 
two meson loop function containing the two pseudoscalar 
mesons of the Z — th channel: 
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with P the total incident momentum, which in the cen- 
ter of mass frame is (y/s, 0, 0, 0). The Gi integral above 
is divergent, and hence it has to be regularized, which 
can be done either with a three momentum cutoff, g max , 
or with dimensional regularization. The connection be- 
tween both methods was shown in Refs. @, [ll|. We use 
the cutoff regularization since this was the one used in 
ref. [HI with a value <7 max = 1 GeV which satisfactorily 
describes a wide phenomenology for meson-meson inter- 
action (T^ |. 

In Eq. (JIJ V, the kernel of the BS equation, is a matrix 
containing the s-wave projected scattering amplitude for 
the two pseudoscalar mesons. Since the main source of 
isospin violation is the different mass between the differ- 
ent charge states, specially the kaons, it is much more 
convenient to work in the charge, or particle, basis in- 
stead of the isospin basis used in ref. (34| and in typical 
chiral unitary approach works. Thus the different chan- 



nels considered are K + K , K Q K°, 



r\r\ and 



Tr°rj. 

Like in ref. |33[ the tree level vertices, V, in Eq. (fT]) are 
obtained from the lowest order chiral Lagrangian with 
isospin-breaking/electromagnetic effects [36[, from where 
the potential between the different channels can be evalu- 
ated. The explicit expressions can be found in Eq. (D2) of 
ref. (33j but all the amplitudes have a difference sign with 
respect to our notation and the amplitudes involving a 
7r°7r° or 7777 must be multiplied by an l/y2 factor to agree 
with the unitary normalization of ref. . As explained 
in ref. [33], these elementary meson- meson amplitudes 
express the isospin breaking effects in terms of the Tr°r] 



K+ 



mixing angle e, and the charge to neutral pion and kaon 
mass difference, An = m 2 ,—m 2 and AK — m 
respectively. 

The use of these tree level amplitudes and the different 
masses between different charge states in the two meson 
loop functions is what generates the isospin violation in 
the unitarized amplitudes. 

For an illustration of this effect we show in Fig. [T] 
the modulus of the amplitudes 1—0 < KK\t\int >/=o 
/= i < KK\t\irrj > /=1 and j =0 < KK\t\irq > I=1 . In 
the conserving isospin amplitudes we see clear peaks cor- 
responding to the /o and ao resonances for the isospin 
1 = and 1=1 respectively which are dynamically gen- 
erated. The i = o < KK\t\irr] >i=i amplitude manifests 
the fo — ao mixing dynamically obtained within the chi- 
ral unitary approach. Note the difference between the 
strength of the isospin violating amplitude with respect 
to the conserving ones which gives an idea of the size of 
the isospin violation through fo — ao mixing. 
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FIG. 1. Modulus of the isospin conserving amplitudes 
i=o (KK\t\wir) / _ Q , /=i (KK\t\w7]) / _ 1 and the isospin violat- 
ing i =0 (KK\t\KT}) I=v 



III. THE J/* -» 07r°77 MODEL WITHIN THE 
CHIRAL UNITARY APPROACH 



In this section we explain the model we use for the 
evaluation of the J/^ — ► <fiir°r), which is based in our pre- 
vious work of ref. [34| , and how to implement into it the 
isospin violation explained in the previous section. Our 
model of ref. [34[ successfully identified and established 
the relevant mechanisms that contribute to the decay of 
the J/^f into one vector meson and two pseudoscalars. 
In particular, ref. [34] focused on the isospin conserv- 
ing channels ujtttt, c^tttt, uKK and <j)KK. The main 
contribution was found to be the mechanisms contain- 
ing the "direct" J/^VPP vertex implementing also the 
final state interaction of the pseudoscalar pair using the 
techniques of the chiral unitary approach from where the 
scalar mesons appear dynamically in the meson-meson 
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scattering amplitudes without the need to include them 
as an explicit degree of freedom. The actual shape of 
the amplitudes in the complex plane, and in particular 
in the physical real axis, couplings of the scalars to the 
different channels, etc, come naturally with the only in- 
put of the lowest order chiral Lagrangian and the im- 
plementation of unitarity in coupled channels. Only one 
free parameter (for the regularization of the two meson 
loop function) is fitted to the experimental meson-meson 
data. (See for instance refs. jg, HU and the references 
cited in the introduction). The model also includes other 
mechanisms where the J/^S decays into a vector and or 
axial-vector and a pseudoscalar meson and the vector 
subsequently decays into the final vector an the other 
final pseudoscalar. Also the final state interaction is im- 
plemented as before. These new mechanism where found 
to be crucial in order to obtain the actual shape and 
strength of the mass distributions in the different chan- 
nels because, in spite that they are small by themselves 
the interference with the dominant mechanism is very 
important. 

In this section we extend and adapt the model of 
ref. [13] to the J/ty — > (f>ir°7] stressing the differences of 
this channel with those studied in ref. [34|. We refer the 
reader to [13] for the details and here we only address 
the particularities, differences and the points of special 
interest for the present work. 



processes involving isospin violation. We do not con- 
sider tree level Jf^^n ^ vertex since in ref. [H| was 
shown that isospin violation in the production operator 
is very small, of the order of 4%. Other possible source of 
isospin violation like the soft photon exchange in the me- 
son propagation was also found to be small in ref. [H[. 
Thus the main contribution to isospin violation comes 
from the /o — ao mixing which in our formalism is gener- 
ated essentially from the loops and potentials implicit in 
the meson-meson unitarized scattering amplitudes, thick 
dots of Figs. H and H 

For the evaluation of the diagrams of Fig. [3J which we 
call "direct" production in the following, we need first 
the J/^<pPP vertex which in ref. [H| it was constructed 
using SU(3) arguments to relate the different channels 
and also parametrizing the amplitudes in a way which 
manifest the OZI rule violation, in a similar way as in 
ref. [13] ■ This introduces two a priori unknown constants: 
an overall coupling g and the OZI violation parameter, 
Xtf,. These free parameters were fitted in ref. |34j to the 
invariant mass distributions of the J/ty — > unt + Tt~ and 

Following the procedure of ref. [34| we obtain 
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After the implementation of the rescattering of the fi- 
nal pseudoscalar pair, the amplitude for the direct pro- 
duction is given bjQ 



FIG. 2. Diagrams for the final state interaction from the di- 
rect J/tycpPP vertex implementing the final state interaction 
represented by the solid circle. 
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FIG. 3. Diagrams for the sequential exchange of vector 
(if* (892)) and axial-vectors (ifi(I270) and ifi(I400)) with 
final state interaction. 



The J/\E r —> (/m°r) decay can proceed only through 



The mechanisms with sequential exchange of vectors 
and axial-vector mesons for the J/^S — > <fiir°r] decay 
are shown in Fig. [3] For this particular decay the 
vector exchanged is a if* (892) and the axial- vectors 
are the K i (1270) and K i (1400). We need the vertices 
J/'i'-vector-pseudoscalar (J/^fVP), vector- vector- 
pseudoscalar (VVP), J/^P-axial-vector-pseudoscalar 
(J/^BAP) and axial-vector-vector-pseudoscalar (AVP) 



1 In Eq. J3J the amplitudes involving 7r°7r° or rfq are in good nor- 
malization, not in the unitary normalization used in the Bethe- 
Salpctcr equation, Eq. Q. 
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which we obtain from the Lagrangians of refs. [34|, |38| . 
[U and respectively. After implementing the final 
state interaction of the pseudoscalar pair the amplitudes 
for the sequential exchange of the K* meson is 



tr = GGt 1 ^ (t K+K - 



(4) 



where G is the coupling in the VVP Lagrangian, G the 
coupling in the J/^VP one, and contains the three 
meson loop function, K* KK, including the momentum 
structure of the different vertices. (See ref. 34] for details 
and values of the constants^) . 

For the ifi(1270) axial- vector exchange we get 



M J/^ m ^ m Kl (1270) 



cD{cD + sF) 



X {tK+K- 



(5) 



,7T°J) + t K°K°,TT a r,) 

where D and F are couplings in the AVP Lagrangian, 
D is the coupling of the J/^AP, s and c parameterize 
the mixing between the isospin 1/2 members of the axial- 
vector octets to give the physical if i (1270) and if i (1400) 
states and i'^ v contains the three meson loop function, 
ifi(1270)ifif , including the momentum structure of the 
different vertices jjyj]. For the diagrams with ifi(1400) 
intermediate state the amplitude is the same but chang- 
ing m Kl(1270) -> m Kl(1400) , F -> -F, c->sands->c 
and replacing the masses and widths of the if i (1270) by 
those of the If i (1400) in the evaluation of i'"" '. 

IV. RESULTS 

The main observable that we are going to evaluate is 
the differential decay width of the J /iff with respect to 
the ir°r] invariant mass, M n o„, which in the J/ty rest 
frame can be evaluated as 



To obtain the total scattering amplitude, t = e*e v t^ v , all 
the mechanisms explained in the previous section must 
be added. 
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FIG. 4. Different contributions to the Tv°ri mass distribution 
in the J/\E' — > (f>-K ri decay. In the theoretical calculation 
we have added the experimental background, dotted line, as 
given in [3lJ to ease the comparison with the data. 
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where is the step function, cos( 



(6) 



^») 2 - m 2 , - |<f | 2 - |Aro| 2 )/(2|<f (/ ,||p 7r o|), where 9 is the 

angle between p* n o and q$ and u>i — \/qi 2 + mf 
hnal polarization sum and initial average is 



The 



<i 2 = 5 E 



(7) 



2 In ref. [341 ] two different results for the fit were obtained depend- 
ing on the sign of G. Our result in [33l for G > is compatible 
with the later work of ref. [4(j | and thus we use this result of the 
fit of |34| for the evaluation of our central results. Our other 
result with G < will be included in the error uncertainties. 



FIG. 5. Same as Fig. [4] but zoomed around the peak region 

In figs. |4] and [5] we show the contribution of the dif- 
ferent mechanisms to the ir°7i mass distribution, -j^- — . 

(Fig. [5] is just a zoom of Fig. 2] around the peak region). 



The experimental data is the result from BESIII [31| 
The dotted line represents the background given in the 
experiment [3l| . Since it is not subtracted in the origi- 
nal experimental data, we have added this background 
to our theoretical results shown in the figure for the 
sake of comparison with the BESIII data. The data 
provided in ref. [3l| is given in number of events but, 
since ref. [31( also provides the total number of J/ ^ pro- 
duced and detector efficiencies, it is possible to obtain 
the absolute normalization for the observable shown in 
the figure. The solid line represents our full theoretical 
result. The dashed line represent the contribution of the 
direct J/^VPP mechanism, fig. [5] The dashed-dotted 



5 




950 960 



970 980 990 1000 

M o (MeV) 



FIG. 6. Solid line and error band: same as in fig. [5] Dashed 
line: full result plus experimental background implementing a 
folding with the experimental resolution, (bins 10 MeV wide). 



line is the result from the mechanisms with sequential 
exchange of the vector meson, upper row in fig. G2 and 
the dashed-double-dotted line is the result from the se- 
quential exchange of axial- vectors, bottom row in fig. [3] 
We see that the loops from direct mechanisms by them- 
selves, which is essentially the only mechanism considered 
in ref. [33| overestimate the strength in about a factor two 
with respect to the experimental data at the peak. The 
new mechanisms containing the sequential exchange of 
the vector meson (K*) and the axial- vectors, (ifi(1270) 
and l^i (1400)), are not very large by themselves but they 
are crucial in order to provide the right strength of the 
distribution when added coherently to the direct mecha- 
nisms. Thus these sequential mechanisms are important 
in order to obtain the right /o — clq mixing. 

In Fig. [5] we show our final result including the theo- 
retical error band obtained implementing a Monte Carlo 
gaussian sampling of the parameters used in the model 
within their error bounds. (See ref. [34[ for the val- 
ues and uncertainties of these parameters). To ease the 
comparison with the experimental data we also show 
with dashed lines our full result averaged in bins of 
10 MeV, like the experimental resolution. Note that 
there is no fit to this observable. All the parameters 
of the model are obtained from other reactions or the- 
oretical works, thus our results are genuine predictions. 
For the branching ratio BR{J/^ — > r/w ^) we obtain 
(4.8 ±0.5) x 10~ 6 which can be compared to the BESIII 
result (5.0 ± 2.7(stat) ± 1.7(sys)) x 10" 6 [3l|. 

In ref. [4l[ an observable called "mixing intensity", 
£/ a , was defined in order to quantify the effect of the 
/o(980) — > a (980) transition independently of the pro- 
duction reaction. For the present case it would be defined 
as 



dl 



1/ * -> tpfo ->0O O 7] 



dr 



J/*-><fi/0-K7i7r71- 



(8) 



in the final state. In our formalism it has no sense to 
strictly isolate the /o or ao contributions by themselves 
since the chiral unitary approach generates the full am- 
plitude from where the definition of the resonance is a 
matter of convention. Thus, the meaningful observable 
is the ratio 



J/*->07r ?; 



J/'Sf—KpTTTT 



(9) 



which close to the /o mass region should be very similar 
to Eq. ([5]). The calculation for the J/^f —> (jmn within 
our formalism was done in Ref. [340. 

The result for the mixing intensity is shown in Fig. [7] 
The mixing intensity can also be obtained from other pro- 
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FIG. 7. fo — > ao mixing intensity. 



duction processes, like e.g. radiative <j> decays 4> ~ * Tr°rjj 
and (p — > ir°r]j. In order to compare the mixing in- 
tensity with experimental data, in ref. (4l| the mixing 
intensity at M = 991.4 MeV, which is the two kaon 
threshold for average kaon masses, was obtained for 
different experiments (42Tj47| . From Fig. [7] we obtain 
£ /o (M = 991.4 MeV) = 0.020± ;g?f to be compared with 
the experimental results, quoted in ED, 0.088, 0.034, 
0.019, 0.027 for SND gl gj, KLOE [H gj], BNL [H 
and CB [U respectively. 



V. SUMMARY 

We have studied theoretically the isospin violating 
J/^> — > (fin rj extending our previous model [34| on the 
J/^> —> (j>PP decays. The main contribution of the model 
are the mechanisms containing the "direct" J/^VPP 
vertex implementing also the final state interaction of 
the pseudoscalar pair. This final state interaction is 



In Ref. 34] the J/ty — > (fnr + ir channel is provided. One has to 



where M is the invariant mass of the two pseudoscalars 



take into account that T 

gO to 7T + 7T 



and 1/3 to 7t 7t°. 



jTj/*^^^^ since 2/3 
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evaluated using the techniques of the chiral unitary ap- 
proach, which implements unitarity in coupling channels 
with the only input of the lowest order pseudoscalar- 
pseudoscalar chiral Lagrangian. At the relevant energies 
of the present work, the scalar / (980) and /o(980) ap- 
pear naturally (dynamically) in the unitarized scattering 
amplitudes without the need to include them as explicit 
degrees of freedom. Due to the strong interference with 
the previous mechanisms, we find also very relevant for 
the final results the mechanisms where the J/^> decays 
into a vector and or axial- vector and a pseudoscalar me- 
son and the vector subsequently decays into the final vec- 
tor and the other final pseudoscalar, where the final state 
interaction is also implemented. 

The isospin violation comes essentially from the differ- 
ence between the masses of the charged and neutral kaons 
and pions in the loops and the meson-meson scattering 
amplitudes. 

We find a very good agreement with the experimental 
BESIII results for the ir°r] mass distribution and branch- 



ing ratio. It is worth noting that without the inclusion 
of the sequential exchange mechanisms we would have 
obtained results a factor about 2 larger. 

The fo — ao mixing intensity, as defined in ref. fill ], is 
also evaluated from the ratio between the J/$ — > (f>ir i] 
and J/\& — > 07T7T decays and we obtain a fair agree- 
ment with most of the experimental results evaluated in 
ref. [Hj]. 
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